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ABSTRACT 

Recently, new observations of ^Li in Pop II stars of the galactic halo have shown a surprisingly high abun- 
dance of this isotope, about a thousand times higher than its predicted primordial value. In previous papers, a 
cosmological model for the cosmic ray-induced production of this isotope in the IGM has been developed to 
explain the observed abundance at low metallicity. In this paper, given this constraint on the ^Li, we calculate 
the non-thermal evolution with redshift of D, Be, and B in the IGM. In addition to cosmological cosmic ray in- 
teractions in the IGM, we include additional processes driven by SN explosions: neutrino spallation and a low 
energy component in the structures ejected by outflows to the IGM. We take into account CNO CRs impinging 
on the intergalactic gas. Although subdominant in the galactic disk, this process is shown to produce the bulk 
of Be and B in the IGM, due to the differential metal enrichment between structures (where CRs originate) and 
the IGM. We also consider the resulting extragalactic gamma-ray background which we find to be well below 
existing data. The computation is performed in the framework of hierarchical structure formation considering 
several star formation histories including Pop III stars. We find that D production is negligible and that a po- 
tentially detectable Be and B plateau is produced by these processes at the time of the formation of the Galaxy 
(z-3). 

Subject headings: Cosmology - Cosmic rays - Big Bang Nucleosynthesis - IGM 



1. INTRODUCTION 

Big Bang Nucleosynthesis (BBN) together with Galactic 
Cosmic-Ray Nucleosynthesis (GCRN) has revealed a consis- 
tent picture for the origin and evolution of deuterium, he- 
lium, lithium, beryllium and boron. This combination in- 
volves very different aspects of nucleosynthesis including 
primordial, non-thermal and stellar nucleosynthesis, all of 
which are correlated through cosmic and chemical evolution. 
There is one single free parameter in the standard model of 
BBN, the baryon density, and that has now been determined 
with high precision by WMAP ( Spergel et al. 200 ^ rendering 
BBN a parameter- free theory dCyburt. Fields. & Ohvdl2002L 
l2003h . The derived BBN value of D/H is in agreement with 
that deduced from observations (D/H) = 2.84 ±0.26 x 10"^ 
jO'Meara et al. 2()06, and references therein). This is a key 
success of Big Bang cosmology. 

Unlike deuterium, which is observed in high redshift quasar 
absorption systems, the LiBeB abundances are primarily 
determined from observations of the atmospheres of stars 
in the halo of our Galaxy. Low metallicity stars (Pop II 
stars) offer various constraints on the early evolution of 
those light elements, assuming that time and metallicity 
are correlated. Until recently, the evolution of the abun- 
dances of ''Li, ^Li, ^Be, and "^^''B could be explained in 
the context of GCRN along with the primo rdial value of 
^Li from BBN which forms the Spite plateau (ISpite & Spitd 
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|I982|) for stars with metallicity lowe r than about [Fe/H] 
= -1.5 (for a review see Vangioni-Fl am et alj l20b0i) . The 
first observations of ''Li a t [Fe/H] ~ -2 jSmith et al.llI993l: 



Hobbs & ThorburnI 1 1 994 Il997t ISrnith. Lambert. & NissenI 



1998t ICayrel et alJ 1 1 9991: iNissen et al.1 l2000l) were entirely 
consistent with the predi cted abundances of ^Li/^Li c± 0.05 in 
standard GCRN models (ISteigman et alJl993HFrelds & Olivd 
[1991 1 Vangioni-Flam et al.llT9m 

Unfortunately, recent observations have led to two dis- 
tinct Li problems. First, given the baryon density inferred 
from WMAP, the BB N predicted values of ^ Li are ^Li/H 
= 4.27+^:02 X 10-1" ■ ■ ■ ■ 



ICvburtll200l . ^Li/H 



or 



^Li/H 



JCvburt. Fields. &"oiiva IMTl llool 
4.9^1: ^ X 10-1° (ICuocoet al.ll200l . 
4.15+j] j^ X 10-'" (ICoc etal.ll2004l) . These val- 
ues are all significantly larger than most determinations of 
the lithium (^Li H-^Li) abundance in Pop II stars which 
are in the range I -2 x 10-'° (see e^ . ISpite & SpitdlI982l: 
iRvan et al.l2000l : lBonifacio et alj20"07l). These values are also 
larger than a recent determination by iMelendez & RamfrezI 
(l2004l) based on a higher temperature scale. This is still 
an open question which will not be addressed in this pa- 
per Se cond, recent observations of ''Li at low metallicity 
dAsplun d et al. 2006; Inoue et al. 2005) indicate a value of 
[^Li] = log^Li/H+ 12 = 0.8 which appears to be indepen- 
dent of metallicity in sharp contrast to what is expected from 
GCRN models, and is more consistent with a pre-galactic ori- 
gin. While BBN does produce a primordial abundance of 
^Li, it is at the le vel of about 1000 times below these re- 
cent observations jThomas. Schramm. Olive. & FieldsllI993l : 
IVangioni-Flam et alj|1999l) . GCRN builds on the BBN value 
yielding a ^Li abundance which would be proportional to 
metallicity. Thus, an additional source of ^Li is required at 
low metallicity ([Fe/H] < -2) and different scenarios have 
been discussed which include the p roduction of ^Li dur- 
ing the epoch of structure formation jSuzuki & Inoudl2002l : 



iNakamura et alj2006l:lT^tischefF & Thibaudl2007l ) or through 
the decay of relic part icles during the epoch of the big bang 
nucleosynthesis (e.g. IKaw asaki et alJ 12005 : iJedamzik et alj 
20051: K usakab eJCaiino & Mathewsf Hooet IPospelov 1 120061: 



Cyburt et al.ll2006ir 



In previous papers (iRoUinde. VangiortT & OUvd I2OOI 
I2OO6I hereafter RVOI, RVOII), we investigated the possibil- 
ity of high-redshift Cosmological Cosmic-Rays (CCRs), ac- 
celerated by the winds of Pop III SN, and thus related the 
non-thermal production of ^Li to an early population of mas- 
sive stars. The star formation history was computed in the 
framework of the hierarchical stru cture formation scenario, 
as described in lDaigne et alj (l2006h . RVOII showed that the 
production of ^Li by the interaction of CCRs with the IGM 
provides a simple way to explain the observed ^Li abundance, 
within a global hierarchical structure formation scenario that 
accounts for reionization, the star formation rate (SFR) at red- 
shift z < 6, the observed chemical abundances in damped Ly- 
man alpha absorbers and in the intergalactic medium. The ad- 
ditional amount of ^Li produced is negligible, so that the dis- 
crepancy between the theoretical prediction and the observed 
Spite plateau is not worsened. 

Here, we examine the consequences of the CCR produc- 
tion of ''Li on the abundances of the related Be and B iso- 
topes. The CCR spallation of p, a on CNO in the IGM as 
well as the reverse process of the spallation of CNO in CCRs 
on H and He in the IGM are considered in computing the 
IGM abundances of the LiBeB elements. As a consistency 
check, we also compute the abundance of D in the IGM and 
the extra-galacti c 7-ray background (EGRB) produced by the 
same processes (Silk & Schramm|[l992l: [Fields & Prodanovid 



l255llPavlidou & Fieldsll2002ir As in RVOII, we assume that 
the IGM abundances act as a prompt initial enrichment (PIE) 
for the halo stars when the galaxy forms (at z ^ 3). The 
observed ^Li plateau at low metallicity is assumed to origi- 
nate solely from production in the IGM by those cosmologi- 
cal CRs, which sets the single free parameter, taken to be the 
fraction of energy available for CCRs. 
In addition, we include two low en- 

lOlive et all 



ergy com ponents, as proposed by 



(119951) : 

process 



(Il994l): ICasse. Lehoucq & Vangioni-FlamI 

IVangioni-Flam et aP (1 19961): t he n eutrino ^ 

dWooslev et al.l 119901: iHeger et al.] l2005h and low-energy 
interactions of C,0 and a with the ISM (LEC). Both should 
be operative in high redshift structures and produce LiBeB. 
These processes are, however, confined within the environ- 
ment of the star and affect mainly the ISM, while the standard 
CRs escape efficiently into the IGM. 

In Section |2] an overview of the basic CCR s cenar io is 
given. We describe the IGM production of LiBeB (§ 12. Il l, dis- 
cuss the escape efficiency of CCRs (§ 12. 2b . and our inclusions 
of the primary z^-process and LEC components (§ 12.31 1. In 
Secti on [3 ] our res ults f or the IGM abundances of ^Li (§ 13. 11 1. 
D (§ 13. 2b . BeB (§ 13. 3b and 7-rays (§ 13.4b are given, and then 
summarized in Section|4] The assumed cosmology is ACDM 
defined by 17,,, = 0.3, 17a = 0.7, h = 0.71 and Q.hlr = 0.0224. 

2. THE CCR SCENARIO AND PRODUCTION MECHANISMS 

Our modeling of the nucleosynthesis by CCRs is based on 
the global scenario for hierarchical structure growth an d cos- 
mic star formation, as developed by Da igne et al.l (l2006f). This 
model is based on an Press & Schechter formalism and satis- 
fies a number of observational constraints such as star for- 



mation rates at z < 6, reionization at high z, as well as the 
abundances of several trace elements in the ISM and IGM. 
Once the initial mass function (IMF) and SFR are specified, 
one can compute the supernovae driven outflows which enrich 
the IGM. The same supernova rate will be used here to derive 
the flux of CCRs. The efficie ncy f or the escape of these CCRs 

is discussed below in Sectio nl2.2l 

The models considered in lDaigne et alj (|2006|) were all bi- 
modal models of star formation. Each model contained a nor- 
mal mode of stars with masses between 0.1 Mp, and 100 M© 
and an IMF with a near Salpeter slope. The SFR of the nor- 
mal mode peaks at z w 3. In addition to a normal mode of 
star formation there is a massive component which dominates 
star formation at high redshift. Here, we consider two possi- 
bilities for the massive mode. First, as in RVOII, we consider 
a model where the massive mode corresponds to stars with 
masses in the range 40-100 M0 (so-called model 1). These 
stars terminate as type II supernovae. Second, we consider 
a model where the massive mode corresponds to stars with 
masses in the range 270-500 Mp (so-called model 2b). These 
massive stars are assumed to terminate as black holes through 
total collapse and do not contribute to any metal enrichment in 
either the ISM or IGM. It is, however, unclear whether these 
implosions are responsible for the acceleration of cosmic rays. 
Energy must get out during the collapse, but this may be en- 
tirely in the form of neutrinos and gravitational waves. We 
have not included any contribution to the flux of CCRs from 
the massive component of model 2b stars. 

2.1. CCR origin 

CCRs are made by elements initially present in the ISM, ac- 
celerated by the SN shock-waves, and then expelled from the 
structures. Nucleosynthesis occurs when they finally interact 
with the elements present in the IGM. 

We compute the energy and flux of CCRs as in RVOI and 
RVOII. The fraction of the total energy injected in CCRs is 
ccr/IOO (where 99% of the energy is emitted in the form of 
neutrinos), with an efficiency ecR for CR acceleration. The 
total kinetic energy S' per SN in CRs is, Sji = lO^^'^ecR ergs 
for stars which leave neutron stars as remnants, i.e. stars with 
masses, 8 M© < m < 30 Mg, and are associated with Pop II. 
For more massive stars, up to 100 M,;,, we take the total en- 
ergy of core collapse to be 0.3 times the mass of the He core, 
with Mhc = 5^ • (OT-2OM0) (IHeger et al.l2003l) . Thus, for the 
massive mode (associated with pop III) of model 1, which is 
dominated by 40 M.j stars, we have — 10^^'^ecR ergs. In 
RVOII, ecR, which within reason is a free parameter for each 
model, was taken to be 0. 15 for model 1 . 

2.1.1. CCR Spectra 

As in RVOI and RVOII, the source spectra are taken to 
be power laws in rnomer itur n, dQ/dp °^ p~'^ dDrurvl [T983I : 
iBlandford & Eichleill 19871 or lJonesli 19941 for a brief introduc- 
tion). The CR proton flux is normalized to the total kinetic 
energy injected by the SN. 

In our scenario, we have considered a simplified picture for 
CCR generation, where the source and propagated spectra are 
equal. As in the RVOI and RVOII analysis, where the CR 
spectral index 7 = 3 was assumed this is consistent with the 
limited scope of our scenario, where the details of injection, 
propagation in the structure, escape, and propagation in the 
IGM are advantageously described by a single phenomeno- 
logical source/propagated power law spectrum. We will come 
back to and comment further on this issue in our conclusions. 
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Given this phenomenological approach, we choose to stay 
as close as possible to existing data. Because of the lack of 
a sound description of CR spectra at all redshifts, a possi- 
ble and conservative choice for describing the CR spectra of 
heavy elements (CNO) is to match the observed present-day 
Galactic CR spectrum. Indeed, a departure from a pure power 
law is observed at low energy: this is caused by the well- 
understood galactic confinement and preferential destruction 
of heavier elements co mpared to lighter ones (see, e.g.. Fig. 1 
of iMaurin et al.l |2004|) . For simplicity, we assume a pure 
power law for p and a with spectral index 7p and 7He- In ad- 
dition, the spectrum for CNO is modeled as a broken power- 
law: 



dp 



■ p-lpM. if £ > £o 

p-7c,N.o otherwise . 



(1) 



Based on observations of the H, He and CNO spectra (see 
Fig. 13 in iHorandei l2007h . we set 7p ?» 7He = 7 3, and 
7c ~ 7n ~ 7o ^ 1 -5 with £0 = 8 Ge V/n. Such a parametriza- 
tion will prove to be crucial to the resulting BeB production 

(§[33). 

In the following, all CCRs above a given energy Ecut = 
1 Me V/n are assumed to efficiently escape from the structures. 

2.1.2. CCR and IGM abundances 

The abundance pattern in the CRs is linked to the abun- 
dance in the ISM. In Galactic Cosmic Rays (GCRs), signif- 
icant deviations for some elements are observed due to dif- 
fering acceleration efficiencie s (e.g. lEllison. Drurv & Mevej 
[1997; Meyer, Drur y & Elliso ni ll997h . The same behavior is 
assumed for CCRs. It is convenient to define the abundances 
of different species /={He, C, N, O} relative to H in CCRs, 
the ISM and the IGM: 



f _ "iSM 
// — H 

"ism 



and 



_ "IGM 



'IGM 



(2) 



where <&, is the CCR flux for species / and n' is the respective 
number density. The quantity F, is related to the interstellar 
abundances by, 

Fi^Vifi, (3) 

and rji is matched to the abundance pattern observed 
in the Galactic Cosmic R ay fluxes (see Table 1 of 
iMever. Drurv & Elhsonl 19981) : 7?He = 0.7, 77c = 8.45, 770 = 5 
and ?7N = 1.47. Abun dances in the I SM (f j) and in the IGM 
(fi) are computed as in lDaigne^t al.l (120061) . In the following, 
we will always assume /hc = = 0.08. The abundances of 
metals, on the other hand, /c,n.o and /c.n.o, do depend on z- 
Both Fi and ft are ingredients to the IGM production and 
their evolution with redshift for C and O are displayed in 
Fig-fflfor the two SFR models considered. Note that the con- 
tributions of reactions involving the CCR flux of N (not shown 
in the figure) are always negligible compared to C and O. Note 
that the CR fluxes and IGM enrichment in model 2b evolve 
more slowly than in model 1, due to the presumption that 
the very massive Pop III stars associated with model 2b, pro- 
vide neither metal enrichment nor a source for cosmic rays. 
In model 2b, both originate from the normal mode, whereas 
in model 1 , CRs and metal enrichment receive contributions 
from the normal mode and massive mode at higher redshift. 



10+0 



u 



10- 



e 
s 



10 



-10 



: Model 1 
: Model 2b 




10 



20 



30 



redshift,z 



Fig. 1. — Relative abundances of C and O in CCRs (f): colored, thick 
lines) and in the IGM (J-,: black, thin lines) for model 2b (solid) and model 1 
(dashed). F; and f; are defined in Eqs. j2) and (3)- 



2.2. Cosmic Ray escape efficiency 

In RVOI and RVOII, the efficiency of secondary CR escape 
from the ISM of early galactic structures (hereafter simply 
"galaxies") into the IGM was always assumed to be high. 
However, as shown in RVOII, the effects of CR heating of the 
IGM, required a certain degree of confinement of CR prop- 
agation and ^Li production to regions that will develop into 
the warm-hot IGM. Here we present a physical justification 
of these assumptions. 

In accordance with iDaigne et al.l ( 120061) . we can evaluate 
the mass M{z) of a typical galaxy at redshift z as the mass- 
weighted average over the Press-Schechter mass function of 
collapsed dark matter halos fps{M,z), 



dMM^fps{M,z)/ / dMMfps{M,z), (4) 



M{z) 



where Mn,in = lO^M©. The average ISM density and radius 
of the galaxy can be estimated respectively as 

(5) 



and 



R{z) = (3M(z)/47rp,(z)A,(z)) 



1/3 



(6) 

where pc{z) is the critical density of the universe at z and 
At.(z) is the density contrast of halos virializing at z (e.g. 
Barkana&Loeb2001). 

CR escape out of the ISM can be mediated either by dif- 
fusion or by advection in a galactic outflow. The CR diffu- 
sion coefficient n is governed by magnetic fields within early 
galaxies and is quite unce rtain. Here we follow the plausi- 
ble physical prescription of lJubelgas et al.l (120061) and assume 
k{p,z) = 3 X lO^^cm's"' (p/mpc)'''^(«isM/cm-^)-'/^, where 
p is the CR momentum. Note that we have normalized to the 
present-day galactic value at unit ISM number density nisM, 



4 



and neglected the weak dependence on ISM temperature. The 
timescale for diffusive escape is 

Tw{p.z)=R{zf/^n{p,z). (7) 

Alternatively, CRs may escape by being advected in 
supernova-driven outflows of ISM gas into the IGM. The 
timescale for advective escape is simply 

Tadv(z)=/^(Z)/V esc(z), (8) 

where Vesc(z) can be evaluated as in lDaigne et al.l (l2006l) . 

Ionization losses within the ISM may deplete LiBeB- 
producing, low energy CRs before they manage to escape. 
The ionization loss timescale for a CR particle of velocity /3c, 
charge Z and mass Anip in a medium of neutral H is 

Tion(/3,z) = (Ampm,cV4^z2/)/?(7- 1) (9) 

x[ln(2m,cV/3V/H)-/3']"'«rsM, 

where 7 = (1-/3^)"'/^ is the particle Lor entz factor and /h = 
13.6 eV is the H ionization potential (e.g. lLongairlll992h . 
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Fig. 2. — Timescales for CR escape from typical galaxies into the IGM 
due to advection r^jv (dotted, blue) and diffusion rjjff (dashed, red), and for 
ionization losses within the ISM tIo,, (dot-dashed, green), for a particles of 
energy 30 MeV/nucleon, compared with the age of the universe T^je (solid, 
black). 

In Fig. [2] Taijv, Tdiff and rjon, the latter two for a particles 
of kinetic energy 30 MeV/nucleon, as functions of z are com- 
pared with the age of the universe Tage, a rough measure of 
the available time at z- We can see that the timescales for CR 
escape out of the ISM by either diffusion or advection are al- 
ways much shorter than Tage at all z, with diffusion dominating 
over advection above z^l ■ This justifies the assumption of 
RV02 regarding efficient CR escape into the IGM. However, 
we see that in some cases ionization losses in the ISM can 
become important on the escape timescales, and such effects 
should be included in future, more detailed studies. 

While the timescales for advection and diffusion indicate 
that CRs eventually escape to the IGM, it is important to con- 
sider the spallation timescales to determine whether there are 



any potential effects on the CR spectrum. A simple estimate 
of the spallation timescale, Tjpai ~ [nisM/3ccrtot] shows that 
spallation indeed occurs on timescales shorter than those as- 
sociated with escape. Noting that the total destruction cross- 
section is roughly constant above a few hundred MeV and 
using CTjoj ^ 40 mb, <Ttot ^ 80 mb, and <Ttot ^^300 mb, one 
finds time scales on the order Tspai < 10^ yrs. This implies 
that within structures and certainly within volumes associated 
with the warm-hot IGM, we expect the effects of confinement 
to play a role and hence allows us to assume a spectrum with 
a propagated slope of 7 = -3 for p's and a's and a flattened 
spectrum of the form given in Eq. ([TJ for CNO. 

We note that the above estimate of pisM should become in- 
creasingly inappropriate at lower z when the total halo mass 
approaches Milky Way scales and the gas within them col- 
lapses to a thin disk via efficient radiative cooling. Then the 
ISM will have a much higher density within a small disk scale 
height, and much lower density for the remaining halo vol- 
ume. The CR diffusion time out of the disk should be less 
than Tdiff at low z as evidenced by the known CR escape time 
out of the current galactic disk 10^ yr for GeV CRs), and 
that out of the halo into the IGM should also be less by virtue 
of the low gas density in the halo. 

2.3. Low energy component and neutrino spallation 

In the past, it has been shown that, besides standard CRs, 
there is a need for an additional source of LiBeB in the struc- 
tures for at least two reasons: standard CRs do not repro- 
duce (/) the meteoritic "B/'^B isotopic ratio and (//) the lin- 
ear (rather than quadratic) proportionality between Be (and 
B) and Fe in the halo phase. These two constraints can be 
satisfied taking into accoun t two additional sources of LiBeB 
jVangioni-Flam et al.l2000l) ; (/) the neutrino spallation in the 
He and C shells of SN, which synthesizes ^Li and '^B and 
(//) the break-up of low energy nuclei injected in molecular 
clouds which produces all LiBeB isotopes. 

Those additional processes as well as their nucle- 
osynthesis products take place in the structure dur- 
ing SN explosion. They are computed accord ing to 
Olive et al.l(ll994:ICasse, Lehouca & Vangioni-Flaml (fr995); 
Vangioni-Flam et al.l (11996^ The subsequent enrichment 
of the IGM is due to the global outflow as computed by 
iDaigne et aLl(l2006l) . 

3. RESULTS 

For all elements X, primordial abundances are assumed at 
the initial redshift z = 30 of the model. The production of X 
by CCR interactions in the IGM is integrated down to z = 0. 
In addition, the abundance of each element is affected by out- 
flows: these contributions are shown separately in the figures. 
The calculation for D, Be and B production follows closely 
that of the ^Li, as given, for example, in RVOI. The full cal- 
culation is then compared to a simplified calculation (Eq. IA7I 
and values gathered in Tab. lAll App.lAli. Note that through- 
out the paper, E denotes the kinetic energy per nucleon. 

3.1. Lithium 

The abundance of lithium in the IGM increases with de- 
creasing redshift due to the interaction of CR a's with He at 
rest in the IGM (the cross section peaks at 10 MeV/n). Since 
the total energy in CRs is proportional to ecR. the amount of 
lithium produced depends on both 7 and ecR and of course 
on the assumed model for star formation. Note that lithium 
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is also produced through CNO interactions with H and He at 
rest. These are taken into account in the calculation, but the 
contribution to the total abundance of Li is always subdomi- 
nant in this context 02.1.1I ). 

We have constrained ecR to insure a PIE of ''Li with abun- 
dance ^Li/H = 10~"'^ at z = 3. The results for model 1 and 
model 2b are shown in Fig. [3] (dashed and solid lines respec- 
tively). Model 1 was already discussed in RVOII and it was 
determined there that ecR — 0.15 was necessary to produce 
the requisite amount of ^Li. In model 2b, however, PopIII 
stars terminate as black holes, with no assumed production 
of cosmic rays. Thus, the lithium production is delayed 
and a higher efficiency is required, ecR — 0-5- Such a high 
value is consistent with those found in diffusive shock accel- 



eration models ( Berezhko & Volk 1997; Berezhko & Ellison 
1999 |:lBerezh'ko"& V61k.2000. 2006; Blasi. Gabici & Vannoni 



2005[ lEUison & Cassam-ChenaJ 120051: IXang & Jonesll2006t 



Ellison et al.ll2007h . as confirmed by recen t comparison with 



observations of supernov a remnants (e.g. iBerezhko & VoUd 
1200?; "^ toen et al.ll2005h . 

The additional lithium due t o ou tflows from the inner re- 
gions (as discussed in Section I2.3l l is small (dotted lines in 
Fig. O compared to the ^Li production by CR interactions 
with the IGM. Also shown in Fig. [3] is the evolution of ^Li 
assuming the initial primordial value of ^Li/H = 4.15 x 10"'" 
deter mined by BBN at the WMAP value for the baryon den- 
sity (ICoc et al.ll2004h : as already emphasized in RVOI, the 
additional CCR production of ^Li is negligible compared to 
that produced by BBN. 
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Fig. 3. — Nucleosynthesis of both isotopes of hthium through the interac- 
tions of secondary CCR a particles impinging on He in the IGM including 
the effects of outflow. For a given star formation history (taken here to be 
model 1 and model 2b of Daigne et al. 2006) and a slope of the CCR energy 
spectrum (7 = 3), the acceleration efficiency of the CCRs is set to obtain a 
PIE which equals the observed plateau *Li/H=10""^ at i = 3: this leads to 
^CR = 0.15 for model 1 and ecR = 0.5 for model 2b. The individual contri- 
butions to the lithium abundance due to CR's and outflows from the ISM are 
shown by the thin curves as indicated in the legend. 



3.2. Deuterium 

There are two main production reactions for deuterium, 
namely pH and pHe. For the simplified calculation, assum- 
ing 7 = 3 and the two extreme cases for I{'-f,pi,p2), as given 
in Eq. ( IA6b . we obtain 



D 



pH 



^Li 



: 5 and 



D 



pHe 



^Li 



!50. 



(10) 



This makes a total contribution of about D/^Li ^100 com- 
ing from pH, pHe and the numerically similar reverse process 
ap, independent of the redshift. This is consistent with the 
constant value of 70 derived in the exact calculation. Thus, 
the additional production of deuterium by CCRs at z = 3 is 
about 3.1 X 10"''' and is therefore well below the observed 
abundance D/H~ 2 x 10"^ (O'Me ara et al...2006,) . 

3.3. BeB 

Beryllium and boron are produced via spallation interaction 
between protons or a's in CRs and CNO elements in the IGM. 
The reverse reactions, which correspond to CNO CRs spalling 
on H and He of the IGM, is in fact dominant. This can be un- 
derstood as the metallicity in the ISM (where CRs originate) is 
always larger than that in the IGM. In contrast to the galactic 
disc, the heavy nuclei component of CCRs becomes impor- 
tant due to the metal deficiency of the IGM. Qualitatively, the 
ratio of the forward to reverse processes is given by 



(BeB)po 
(BeB)oH 



fo 
Fo 



[O/H 



IGM 



?/0[0/H]iSM 



(11) 



where /, F and rj are defi ned in Eq. ^ and Eq. (O. From the 
calculated abundances of lDaigneet a l. (2006), the metallicity 
ratio between the IGM and CRs is < 10~^ (see eg. Fig.[T]i, so 
that, the forward process is ^ 10"^ smaller than the reverse. 
For the reverse process, Eq. ( IA7b does not apply in principle, 
since Fi no longer factors out of the z integrand. However, 
in the CRs, Fc < lO""* and Fq < 10"^ (see Fig. [B- An up- 
per limit to the production of ^Be, would then be, for 70 = 3, 
^Beon/^Li ~ 0.3. This compares well with the full calcula- 
tion (upper solid line in Fig. |4|. The leading contribution of 
the reverse process is confirmed in Fig.|4]for the full calcula- 
tion in model 2b (upper solid line compared to dashed line) ; 
the same is also true for model 1 . The same reasoning holds 
for reactions involving C, N. The dominant channel is actually 
O-i-H. Nevertheless, all nuclei are accounted for in the full cal- 
culation below. As one can see in Fig. 21 adopting the spectral 
break, i.e. 70 = 1.5 as described in Section l2.1.1l decreases 
the contribution of the reverse process and as a result lowers 
the Be abundance as shown by the thick solid curve compared 
to the upper thin curve with 70 = 3. 

The results from the full calculation of the BeB abun- 
dances for both models 1 and 2b are given in Figs. (|5])- 
Note once again that BeB production in model 2b is 
delayed due to the lack of metal enrichment in the ISM 
from the massive mode in this model. At z = 3, we find 
''Be/H= 10-'^ '' in model 1 (upper panel) and "^Be/H^ 10"'-^ 
in model 2b (lower panel). This is compatible with the ob- 
served abundance at the l owest metallicity, [ Fe/H1= -3.3, 
''Be/H- 10"'^ in G64-12 dPrimas et alj|2000ah . For boron, 
we find B/H = lO""-^ in model 1 and and B/H = IQ-^—^ 
in model 2b. These values are somewhat below the ob- 
served boron abundance at low metallicity B/H^ 1 0"' ''^ at 
[Fe/H]== -3.0 in BD -13 3442 ( iPrimas et alj|2000bl) though 



6 



pm 



M 
o 



-16 



Model 2b 

: O + H 

:p + 




Fig. 4. — Beryllium production by the interaction of CCRs and the IGM for 
model 2b. The direct reaction p+O (dashed Une) is negligible compared to the 
reverse reaction O+H (soHd thick line). The upper thin soUd curve displays 
the 'Be production as a function of redshift assuming an O spectrum with 70 
= 3 as for p and He at low energy The thin dashed curve assumes 7p,Hc = 3. 
The thick solid curve is based on Eq. (TJ and assumes 70 = 1.5 at low energy. 



B/H is as low as 10"'^ at [Fe/H]= -2.9 in BD -23 3130 
dGarcfa-Lopez et aT]|1998h . Note however, that unlike beryl- 
lium abundances, the boron abundances in particular are 
very sensitive to the assumed effective temperature. The 
abundances quoted abo ve were determined using tempera- 
tures based on the IRFM jAlonso. Arribas. & Martinez-Rogeii 
pH). The boron a bundance in BP -13 3442 would be IQ-'"-^ 
jFields et al.ll2005|) had we ado pted the temperature derived 
bv lMelendez & Ramfrezl(l2004 . Of course Li also scales with 
temperature and at higher effective temperatures, the level of 
the ^Li plateau would be raised forcing one to large CR effi- 
ciencies. 

Our results can be cast in terms of abundance ratios of 
''Be/'^Li and B/^Li 



for model 1 and 



Hi 



: 0.020 



: 0.008 



B 

"Li 
B 



iO.20, 



' 0.09 



(12) 



(13) 



for model 2b. These predicted ratios are independent of the 
efficiency ecR and can be compared to the few examples 
where observational determinations of ^Li and Be exist in 
the same star Restricting our attention to only those stars 
with [Fe/H] < - 2.6, we have Be abundan ce measurements 
in two stars, G64- 1 2 (iPrimas et al.l l2000ah and LP 815-43 
jPrimas et al.l 2000bl) for which there are reliable ^lA PlA de- 
terminations ( lAsplund etall 120061: [inoue et al.ll2005l) . When 



o 
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Modell Model 2b 



: total 
: CRs 
: Outflow 



20 



30 



redshift, z 



Fig. 5. — Nucleosynthesis of beryllium and boron through the interaction 
of secondary CR p's, a's and (CNO) impinging on CNO (H and He) in the 
IGM including the effects of outflow. The evolution of the abundance of 
'Be and " B is displayed as a function of redshift for models 1 and 2b. The 
individual contributions to the beryllium abundance due to CR's and outflows 
from the ISM are shown by the thin curves as indicated in the legend. 



pm 



Model 2b 



o 
h-9 




-14 



: total 
: CRs 
: Outflow 



10 20 
redshift, z 



30 



Fig. 6. — As in Fig.|5] for boron. 



corrected for stellar temperature differences, we find ^Be/^Li 
= 0.005 ± 0.003 for G64-12 and = 0.01 1 ± 0.007 for LP 815- 
43. Unfortunately there are no very low metallicity stars with 
both ^Li and B abundance measurements. 
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The predicted PIE at z = 3 for both elements is of the same 
order as the abundances observed at the lowest metallicity. 
Thus we expect a plateau in both Be and B at low metallicity 
similar to that found for ^Li. Indeed, when the B e abundance 
in G64-12 was first reported (iPrimas et al.l200()al) . the authors 
claimed evidence for a flattening of the observed Be evolu- 
tionary trend at low metallicity. It will be interesting if future 
determinations of Be at low metallicity confirm the existence 
of a Be plateau. 

We have also computed the production of BeB through 
LEC nucleosynthesis and neutrino spallation via outflow to 
the IGM. In all cases, the abundance produced by these pro- 
cesses is less at small redshift, though the abundance of '^B 
produced by the i^-process is not negligible at higher redshift. 
These results are displayed by the thin lines in Figs. (ISll-®. 

3.4. Pionic production of gamma-rays 

The interaction of the CCRs with the IGM also produces 
gamma- rays. Our calculation closely follows that given in 



iPavlidou & Fieldsl (l2002h . Assuming isotropic production, 
the extragalactic differential flux of 7-rays at Earth in the flat 
ACDM cosmology is 



dz 



(14) 



The 7-ray source function qj{z,E') is calculated using the 
simplified analytic formula of Pfrommer & Ensshn (2004): 

qy{z,E')^app(l)p{E')niGMiz)C^~°' 



4m_o 
X — ^ X 



(2Ei^ 



2£; 

m_Q 



with C = 2, 6y ^ 0.14a::^' '' + 0.44 and tTpp 



(15) 



32 X (0.96 + 



A.4_2.4o 



' ) mb. The same input fluxes and IGM densities are 
taken as for the D and LiBeB calculations. Fig.|2]shows that 
the total 7-ray flux produced by CCRs is negligible compared 
to the observed EGRB. 

4. DISCUSSION AND CONCLUSIONS 

Our initial motivation to study CCRs at high redshift was 
to find a physical origin for the anomalously high abundance 
of ^Li observed in metal poor halo stars. A cosmic history 
for structure format ion that reproduces standard observations 
(iDaigne et al.l2006h provides enough energy in SN to produce 
a ^Li plateau at the level of log[^Li/H] w -11.2 via a + a in- 
teractions of CCRs with the IGM. Note that the ^Li produced 
by outflows (via neutrino spallation and LEC) is negligible. 

In this paper, we considered two models of star formation 
histories differing by their Pop III stars. The massive compo- 
nent of model 1 (40- lOOMg) injects a lot of energy in CCRs 
so that the necessary fraction ecR is low 0.15). On the 
other hand, the stars associated with the massive component 
of model 2b (260-500M;j ) are assumed to collapse into black 
holes without injecting, a priori, any energy in CCRs. Con- 
sequently, the fraction ecR associated with SN in the normal 
mode must be higher (~ 0.5). This figure would be reduced if 
super massive stars produce CRs. 

We confirm, as in RVOII, that the ^Li primordial abundance 
always dominates any additional production by CCRs. The 
same conclusion holds for deuterium. Similarly, the observed 
EGBR is much larger than the intensity of photons produced 
by the decay of pions produced by CCR proton collisions. 



EGB Data 



Pionic cont. 




10.0 



100.0 



Ey(GeV) 



Fig . 7 . — Energy spectrum (E^I versus energy in GeV) of photons produced 
by pion decay in th e IGM. Tlie models ar e well below the observed EGBR 
data l lStrong. Moskalenko & Ptuskinl20"07l) . 



In the same cosmological context, we explored the produc- 
tion of other light elements (BeB). To that end, it is necessary 
to track the abundances of metals (CNO) in the ISM and IGM. 
We have shown that the reverse process (i.e. CNO CCRs on 
H and He IGM gas) is the dominant channel to synthesize 
BeB. This is in contrast to the production of BeB in the galac- 
tic disk, where the reverse process co ntributes roughly 20% 
(iMeneguzzi. Audouze & Reeveslll97ll) . This is easily under- 
stood as the metal enrichment in CCRs is inherited from struc- 
ture abundances, which are far higher than those in the IGM. 
Note that we also checked that the neutrino spallation and 
LEC processes are negligible in the IGM BeB budget at the 
time of the galactic formation, i.e. z ^ 3. In all models con- 
sidered, we have shown that BeB synthesized at z ^ 3 is at 
the level of the observed abundances in the lowest metallic- 
ity stars. This is the first theoretical indication of a plateau 
for these elements which does not resort to exotic models of 
BBN. 

Note that these results have been obtained with several as- 
sumptions about the CCR spectra. To some extent, our mod- 
eling, be it for light elements (pHe) or metals (CNO), is lim- 
ited. Indeed, two fundamental ingredients for the calculation 
are i) the low energy form of the source spectra for protons 
and ii) the propagated fluxes to plug in the IGM. The first 
item is crucial for determining the acceleration efficiency ecR- 
For example, changing 7 « 3 to 7 « 2 — hence taking a more 
conventional spectral index for the sources — would lead to an 
unphysical value of ecR > 1 ■ However, we expect that prop- 
agation effects on scales of order the warm-hot IGM would 
lead to a steeper spectrum (i.e. 7 = 3 for p's and a's) and a 
flatting of the CNO spectrum of the form we have assumed 
here. The second item involves several issues that are inti- 
mately connected: the details of CR escape and confinement 
between the structures, the warm-hot IGM and cooler IGM is 
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TABLE A 1 

Simplified description of cross sections used in the approximate calculation 



^Li and ^Be are observed simultaneously. Finally, we note 
that B can be observed directly in high red shift objects, as 
achieved for a damped Lya system at z = 2.6 ( iProchaska et alj 
l2003h . Observations at even higher z may be possible through 
absorption lines in gamma-rays bursts (J. Prochaska, private 
communication), which should provide very valuable con- 
straints on early LiBeB production scenarios. 



related to whether or not the propagated spectrum displays a 
spectral index close to the standard source index (7 « 2) or 
closer to a diffused spectrum (7 w 3) and, more importantly, 
how this evolves with z- Without some degree of confinement 
on the scale of the warm-hot IGM, we would be forced to take 
the same spectrum for p. He and CNO nuclei — as would be 
more natural for a 100% escape to the IGM, and this would 
lead to an overproduction of ^Be. Our simple estimate of the 
spallation timescale within structures indicates that confine- 
ment is indeed playing a role and changes to the CR source 
spectrum will occur 

All of these elements clearly call for a more coherent and 
refined calculation. As just outlined, one issue concerns the 
description of CR propagation in realistic structures, evolv- 
ing with redshift, which also allows for differentiated produc- 
tion in situ and outside structures. Another important issue 
is the possi biUty for heterogeneity in me tallicities at a given 
time (e.g. ISalvadori. Schneider & Ferrara ,2006,) , that can- 
not be handled in t he homogeneous paradigm developed in 
iDaigne et alj ( l2006h . This question could be addressed when 

APPENDIX 

SIMPLIFIED CALCULATION FOR D-BE-B 

The calculation for D, Be and B production follows closely that of the '^Li, as given, for example, in RVOI. In this appendix, we 
estimate the relative efficiency of various channels as well as the relative efficiency of the D, Be and B production with respect to 
the *'Li production. In some cases, the latter is independent of z- Below, E denotes the kinetic energy per nucleon. 

The abundance of the element X relative to H is given by 
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i=CR i=IGM 



:o JE, 



fjafjiEp)F,^piEp,z)dEp 



dz , 



(Al) 



where $p(£'p,z) is the accumulated CR proton flux in the IGM (in proper units), afj{Ep) is the production cross section of X in 
the reaction (/+ j) and /, and Fj are defined in Eqs. (|2]l and Q. Note that ^Li is only produced by CR a's with energy four times 
the final lithium energy. As a consequence, the above formula also applies to [^Li/H], but an extra factor 1 /4 must be added. Two 
simplifications are made: 

• All cross sections <^j{E) are approximated as a constant in the range Ei -E2 and zero elsewhere, i.e. 

cTfj{E)=<jfjXe{E-Ei)eiE2-E). (A2) 
The values of the cross sections for the main processes are given in Tab. lAll 

• $p(£'p,z) is assumed to keep the same power law dependence 7 as the source spectrum at each z, so that it can be rewritten 

$p(£p,z)=A?o(z)$p(£p). (A3) 



as 



If y = H, He, fj and Fi are constants that factor out of the z integrand 

'X' 
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where 



/(7,pi(£i),/?2(£2)) 



-7+1 



-.dp , 



(A5) 



and where p{E) is the CR proton momentum for a given kinetic energy per nucleon E. This function can be integrated as a 
hyperboHc function, but for the purpose of a simpHfied calculation, it is good enough to consider the two limits, that apply either 
when the production is at very low energy (e.g. for ^Li), or at GeV/n energies: 



-7+2 
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We shall now write a simplified expression for the comparison of the production of X with respect to '^Li. The integration over 
z cancels out, and this results in 

-C«;WGM ^^^^ 



"Li 



he iFHe/4) aJi 



where we wrote I{'-f,pi{E[),p2 (£2 ) ) = ^; for short. 
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